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Abstract. For a relativistic system of two scalar particles, we find the Bethe-Salpeter amplitude in 
Minkowski space and use it to compute the electromagnetic form factor. The comparison with Euclidean 
space calculation shows that the Wick rotation in the form factor integral induces errors which increase 
with the momentum transfer Q^. At JLab domain (Q^ = 10 GeV'^ /c^), they are about 30%. Static approx- 
imation results in an additional and more significant error. On the contrary, the form factor calculated in 
light-front dynamics is almost indistinguishable from the Minkowski space one. 

PACS. PACS-key ll.lO.St - PACS-key 13.40.Gp - PACS-key 21.45.-v 



1 Introduction 

The Bethe-Salpeter (BS) equation [T] provides a field- 
theoretical framework for a relativistic treatment of few- 
body systems. It has been extensively studied in the liter- 
ature (see [2] for a review) and used to obtain relativistic 
descriptions of two-body bound and scattering states. 

BS equation is naturally formulated in the momentum 
representation. In Minkowski space, for two spinless par- 
ticles, it reads: 



[(§ +fc)2 



iKM{k,k\p)'pM{k';p) 



(1) 



The interaction kernel Km is given by irreducible Feyn- 
man diagrams. Any finite set of them is an approximation 
of the interaction Lagrangian of the theory under consid- 
eration. Most of works were done in the ladder approxima- 
tion, that is restricting the interaction kernel to its lowest 
order exchange term. Several researches (in particular [3J) 
indicate that the higher order kernels, usually not incor- 
porated in the BS equation, give a significant contribution 
to the two-body binding energy. Studies of these contri- 
butions, step by step, in the BS framework, would be of 
indubitable interest. 

Until very recently, the BS equation had been solved 
only in the Euclidean space, i.e. after performing a Wick 
rotation [J] , in order to remove the singularities due to the 
free propagators. The validity of Wick rotation has been 
proved in [4 for the ladder kernel. For higher order kernels 
{e.g., for the cross box) the possibility of Wick rotation is 
less clear, since one deals with the "partially Euclidean" 



BS amplitude, in which the relative energy fcg = iki is 
imaginary, whereas the total energy po remains real. This 
"partially Euclidean" transformation is a subtle point to 
be checked more carefully. We will see below that it is 
valid also for the cross-ladder kernel. Then the Euclidean 
BS equation (in the rest frame) provides exactly the same 
binding energy as the Minkowski one. 

If we are interested not only in the binding energy but 
also in the electromagnetic (EM) form factors, the Eu- 
clidean BS amplitude in the rest frame is not enough. On 
one hand, when computing the integral for the form factor, 
the rotated contour crosses singularities of the integrand. 
So, the result is not reduced to the naive replacement 
fco = iki which transforms the Minkowski BS amplitude 
into the Euclidean one. On the other hand, form factor 
for non-zero momentum transfer involves the BS ampli- 
tude for non-zero total momentum p. This amplitude can 
be obtained from the rest frame one by a boost, but the 
parameters of this boost for real p and imaginary fco are 
complex. This requires the knowledge of the BS ampli- 
tude in the full complex plane. The continuation of the 
Euclidean amplitude from real axis to the complex plane 
is numerically very unstable and can hardly be done in 
practice. Instead of it, the BS amplitude in the complex 
plane can be found by solving the equation for complex 
arguments. However, the equation in the complex plane is 
no longer Euclidean nor Minkowski one and it is actually 
more complicated than the equation on the real axis. This 
difficulty is avoided in the so called static approximation 
[5], which makes use of the Euclidean BS amplitude only 
but brings an additional error increasing with the momen- 
tum transfer. 

These problems disappear if one expresses the form 
factor through the Minkowski BS amplitude. For the lad- 
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der kernel, the BS equation in Minkowski space was solved 
in . For separable interactions, an approach in Minkowski 
space was developed and applied to the nuclcon-nuclcon 
system in [7]. In [8], the effect of the cross- ladder graphs 
in the BS framework was estimated with the kernel rep- 
resented through a dispersion relation. 

Recently, a new general method to find the Minkowski 
BS amplitude has been developed This method is valid 
for any kernel given by Feynman graphs. In the case of 
spinless particles, it was tested for the ladder and cross 
ladder kernels [TU| . 

Having found the Minkowski BS amplitude, we can cal- 
culate EM form factor without any approximation. This 
allows us to check the validity of Wick rotation in the 
form factor integral, the accuracy of the static approxi- 
mation and, in addition, to make a comparison with light- 
front dynamics (LFD) calculations. This is the aim of the 
present study. A first description of these results, without 
any derivation, can be found in [llj . 

This article is organized as follows. In sect. [2] we briefly 
describe the method [9] for solving the BS equation in 
Minkowski space and give new validity tests. In sect. [3] the 
Minkowski space calculation of form factor is presented. 
In sect. |4] the analogous Euclidean space computation is 
carried out, including its static approximation. In sect. 
[5] the form factor in the LFD framework is calculated. 
The comparison of numerical calculations performed in 
the different approaches considered in this work is given 
in sect. |6l Finally, some concluding remarks are presented 
in sect. [T] 



2 The method 

According to PITU]. the Minkowski space BS amplitude is 
found in terms of the Nakanishi integral representation [21 
[12]: 

■ r-l ^OO 

<pM{k;p) = 7= dz d-y 

V47r J-i Jo 



7 + 



iM2 



k"^ — p ■ k z — iey 



The weight function 3(7, z) itself is not singular, whereas 
the singularities of the BS amplitude are fully reproduced 
by this integral. For example, if we set —-^^gij^z) — 1 
and calculate the integral, we find 

^^^^''■'P^ " [(I + fc)2 _ „i2 + le] [(I - ky - m2 + le] ' 

i.e. just the product of two free propagators in ([T|). <Pm m 
the form ^ is substituted into the BS equation ([T|) and 
after some mathematical transformations one obtains 
the following integral equation for 17(7,2:): 

g{j',z)d'y' 



7' 7 -1- z^im? + (1 — z'^)k? 



di / dz'V{^,z-i,z')g{^\z'), 



where — — jM^ and V is calculated in terms of 
kernel Km [2]. A remarkable point is that equation ([3|) 
is strictly equivalent to the original BS equation ([T]) and 
does not contain any singularity neither in 5(7,2) nor in 
V. Once 51(7, z) is found by solving ([3|), the Minkowski BS 
amplitude is obtained by expression ([2|). 

Equation for the Euclidean BS amplitude is obtained 
from the Minkowski one ([!]) in the center of mass frame 
p = 0, by the Wick rotation and substitution feg = ik^ 



k^ + 7r? 



4 



M^kj 



where 



^!^KE{k,,k;klk')$EiK,k') 



^E{ki,k) = <pM{ih,k;M,p 0) 



(4) 



(3) 



(5) 

and KE{ki,k]k'^,k') = KM{ik4, k;ik'^, k'; M,p = 0). By 
the method developed in [0] we can restore the Euclidean 
BS amplitude, i.e., find the solution of (|4|), by setting in 
(21) fco = ik4 and p = 0. 

All three equations ([1]), ([3]) and ((4|) are equivalent to 
each other in the sense that they give the same M^. In 
[9j the equivalence of ([3|) to the initial BS equation ([!]) 
has been checked numerically for the ladder kernel. It was 
found that the binding energies B = 2m — M coincide 
with high precision. 

For more complete checks, we present here two addi- 
tional tests. In the first one, still for the ladder kernel, 
we compare the Euclidean BS amplitudes obtained by in- 
serting the solution of ([3|) into ([2|), where k = {ik4,k), 
p = (M, 0), with the one found by directly solving ((4|). 
These amplitudes are plotted in fig. [T] (top) . They coincide 
at a level better than 0.2% on the whole range of (|fc|, fc4) 
considered and are indistinguishable from each other in 
the graph. On the contrary, the Euclidean BS amplitude 
strongly differs from the initial Minkowski one, shown at 
the bottom part of fig. [TJ It is worth reminding that so 
different amplitudes (compare top and bottom of fig. [T]) 
correspond to one and the same mass M. 

The second test incorporates, in addition, the cross 
ladder kernel displayed in fig. [21 For a massive exchange 
and a wide range of binding energies B, we carried out 
precise calculations (accuracy better than 0.1%) of the 
corresponding coupling constants a = g'^ / (\Qttw?) both 
by equation ([3]) and by Euclidean equation ([4]) (here g 
is the coupling constant in the interaction Hamiltonian) . 
The results are displayed in the table[Tlin units of m. Their 
coincidence demonstrates the validity of both calculations 
- in Minkowski space by the method 9J and in Euclidean 
one. 

The ladder or ladder -|-cross-ladder kernels are good 
enough to check the applicability of the method [9j. We 
would like to emphasize, however, that both kernels give a 
rather crude approximation of the full interaction. Figure 
[3 shows the ground state mass M obtained, for ^ = 0.15, 
by the BS ladder and ladder -|-cross-ladder kernels to- 
gether with Feynman-Schwinger representation results [3] . 
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Table 1. Coupling constant a for given values of the binding energy B calculated by the eq. ((3} and by the Euclidean BS 
equation (j4| for the ladder +cross-ladder (L+CL) kernel. The exchanged mass: ^ = 0.5. 
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Fig. 1. Top: The Euclidean BS amplitude (for ladder kernel) 
obtained by ([2]) with ko = iki, for different values of k. Or (in- 
distinguishable) the one calculated by direct resolution of BS 
equation in Euclidean space (|4]). The exchange mass is /x = 0.5. 
Bottom: The corresponding amplitude in Minkowski space, ob- 
tained in [9l[T3|. We use the units m = 1. 
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Fig. 3. Ground state mass M obtained by the Bethe-Salpeter 
ladder and cross-ladder kernel, compared with the Feynman- 
Schwinger representation results for an exchanged mass jj, = 
0.15. 

The latter incorporates all the higher order cross box con- 
tributions in the kernel, but not the self energy. Even at 
low binding energies, the ladder and cross-ladder results 
differ by at least 20%, this difference reaching more than 
a factor 2 around B/m = 1. On the other hand, the re- 
sults obtained by Feynman-Schwinger representation de- 
parts strongly from the BS(L-|-CL) as soon as B is bigger 
than 0.05 m. The cross-ladder kernel thus gives a non neg- 
ligible contribution to the total mass of the system in the 
right direction, but a larger contribution remains to be 
included, due to the higher order terms. Notice however 
that the underlying field theory with cubic boson-boson 
interaction is unbounded from below [141 . This instability 
(for the interaction g4>'x^) appears when infinite number 
of the -y^ loops in the field <f) self energy is included [H] . 



p. -p 



3 EM form factor via Minkowski BS 
amplitude 



p,-p +p 



Fig. 2. Feynman cross ladder graph. 



The electromagnetic vertex is shown in Fig. |4l We sup- 
pose that one of the particles is charged. By applying the 
Feynman rules to this graph, we get: 



d^k {p + p' - 2kY 



(27r)4 (fc2 -rrfi + ie) 



4 



J. Carbonell, V.A. Karmanov, M. Mangin-Brinet: Electromagnetic form factor via Bethe-Salpeter amplitude. 




Fig. 4. Feynman diagram for the EM form factor. 



r{^p-k,p)r{^p'~k,p') 

[{p — fc)^ — + if] [{p' — fc)^ — + ie 



(6) 



where r{k,p) is the vertex function, related to the BS 
amplitude by: 



^M{k;p) 



r{k,p) 



[(I + fc)2 _ ^2 + jg] [(I _ -m-^+i€ 



(7) 

Therefore the electromagnetic vertex is expressed in terms 
of the BS amplitude by the formula: 



{p+prFM{Q') = i 



d'^k 

(2^ 



{p+p'-2ky 



X {k^~m^)$M Qp-fc;pj<?A/ Qp'-A:;p'j . (8) 

We multiply both sides of ^ by (p + p')^ and substitute 
in its r.h.s. the BS amplitude in terms of the integral 
So, the form factor is given by: 



(.p+prFM{Q') = 



id^k 



d'ydz d'y' dz' 3(7, z)g{'-/', z') 



2(27r)5 

[{p + pT-'2k-{p+p')] jm^-k^) 
[7 + m2 - iM2 - (ip -kf-p- {\p -k)z- iff 

1 

X ^ 

[7' + to2 - \A'P - {\p' - fc)2 - p' ■ (ip' - fc) ^' - ie] 

To compute this integral, we use the Feynman parametriza- 
tion: 



1 



1 30u2(l - u)2(ilt 



'0 (^au + fo(l-u)) 
and then shift the integration variable: 



fc = fci + i(l + + i(l + z')(l - u) p'. (9) 



Then the integral over d^fci has the form 



• • d^ki 



{kl-c^ief 



(10) 



where c does not depend on ki. Though the calculation 
of this integral by Wick rotation is standard, we explain 
it here in more detail, to emphasize the difference with 
the calculation performed using Euclidean BS amplitude, 
where the Wick rotation cannot be done (see sect. [5] be- 
low). The integrand in pop does not contain linear terms 
in fcio, but only a constant and a quadratic term. It has 
four poles at the values fcio = ±\Jc^ + k\^ie. Their posi- 
tions do not prevent from the counter-clock-wise rotation 
of the integration contour. Therefore, substituting here 
fcio = ik^, we get, for the constant term in kiQ, the follow- 
ing relations: 



(fcf -c + ief 
idk^d'^ki 



dkiod^ki 



(fc2Q - fc2 _c + ie)e 
i2n'^k^dk 2i7r2 



(11) 



[kl + kl + cf Jo (fc2 + c)6 40c4 
Calculating similarly the integral for the quadratic term 
kl d*ki 2i7r2 



{kf -c + ie)^ 



60c3' 



(12) 



we find the following formula, which is exact for a given 
5(7,2): 



1 f°° f 
= 2Wm X 

d-y' / dz'g{j',z') / duu^il-u) 



(13) 



2^-, ^^2 fnum 

f4 

•I den 



with 

fnum = (6C - 5)m2 + [7'(1 -u) + 7u](3e - 2) 

+ 2M2^(1 - + \q^{1 - u)u{l + z){l + z') 

fden =m^+-/'il-u)+JU- M2(1 - 

+ ^Q^{l-u)uil + z)il + z'), 

where = —{p — p'Y > 0- To simplify the formula, we 
use the notation: 

£. = \[l + z)u+^{l + z'){l~u). 



We have also introduced in (|13p the normalization factor 
Nm which is found from the condition Fm{Q) = 1- 



4 EM form factor via Euclidean BS amplitude 

Form factor ([13]) was calculated using a well justified Wick 
rotation in the variable kio, defined by (l9|). As explained 
in sect. O the Euclidean BS amplitude in the rest frame 
<pE{k4; k) is obtained from the Minkowski one (see eq. (O) 
by Wick rotation in the variable /cq- To express the form 
factor through <pE{k4,k), one should make the Wick ro- 
tation, in the variable fco, in integral ^ (for the moment. 
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we ignore the fact that the BS amphtude in ([8]) is not in 
the rest frame; we will come back to this point later). We 
will show that, in contrast to integrals pT|) and p^ . the 
Wick rotation in ([8]) cannot be done without crossing sin- 
gularities. Therefore the form factor cannot be expressed 
through the Euclidean BS amplitude exactly. 

It is enough to illustrate this statement in the simplest 
case, with r(fc,p) = 1 and p' = p = (M, 0), i.e. = 
and V = Q. Integral ([6]) then turns into: 



2i{M - ko) 



That is: 



(27r)4 (fc2 - to2 + ie)[{p - fc)2 - + ie]' 



The first propagator has poles at fco — zty/m'^ + =F ie 
and this does not create any problem, whereas the second 
factor has poles at: 



ko=M± \fn^ 



fc2 =p ■ 



If fc2 < 



both poles are in the r.h.s. half plane 



and the pole at fco — M—^/vi? + k'^ + ie prevents from the 
Wick rotation. The exact result for the form factor should 
incorporate the residue in this pole and therefore it is not 
reduced to the integral obtained from ([5]) by the naive re- 
placement fco = ifc4. If the residue is omitted (or, in realis- 
tic case, if the contributions of other possible singularities 
of r{k,p) crossed by the rotated contour are omitted) the 
result is approximate. In practice, taking into account the 
contributions of these unavoidable singularities is impos- 
sible and, hence, the form factor calculated through the 
Euclidean BS amplitude is always approximate. 

Shifting the variable fco (for example, fco — > —fco + ^Po, 
to transform the argument ^po — ko of the BS amplitude 
in ([8]) into fco) does not help. The situation remains the 
same for non-trivial r{k,p) and for non-zero . 

In addition, there is another reason which does not 
allow to express the form factor via Euclidean BS ampli- 
tude. The latter is determined by eq. ^ in the rest frame 
p = and it is related to the Minkowski one by eq. ([5]). 
However, the form factor is expressed through the BS am- 
plitude with non-zero total momenta p and p' which due 
to scattering are different in initial and final states. Hence, 
after Wick rotation, we need to know 



<PT''{ki,k;p)^$Miih,k;p). 



(14) 



which differs from the Euclidean BS amplitude ^£;(fc4,A;) 
in eq. (O, by non-zero value of p. They are identical only 
at p = 0. The boosted amplitude <P^£°''\ki,k\p) can be 
expressed through <pE{k4, fe), but only for complex values 
of its arguments fc4, |fe|. 

Indeed, the Minkowski amplitude ^Af(fcoi k;p) in r.h.s. 
of (|14p , for real fco and for non-zero p, can be found from 
the rest frame amplitude by a boost. Namely, we can take 
the BS amplitude ^A/(fco, |fc|; M,p — 0) in the rest frame 
and substitute 

ko^ k'o^ j^iPoko -p k), 



Ife'l 



jL/2 _ l2 I jU2 
' • 



^Miko,k;p) = ^Mik'o, Ife'l; M,p = 0). 



To get the Euclidean amplitude, we replace here fcg = 
iki, k'o — ik'^, substitute the result in (|14p and use the 
definition Then the relation has the form: 



#|°°^*(fc4,fe;p) = <Pi5(fc4,fc') 

where 'Psik'^^, k') is the Euclidean BS amplitude in the rest 
frame, depending however on the complex arguments: 



fc/1 



1 



{poki + ip-k), k' 



k\ 



k"^ ~ k' ' 



(15) 



This requires the knowledge of the Euclidean BS ampli- 
tude $E{k'^,k') in the full complex plane. Alternatively, 
one can solve the Euclidean BS equation for non-zero p 
(for real arguments) and obtain ^^°'**(fc4, fe;p) directly. 
These solutions for quark systems were found numerically 
in fTHj. In sect. I4.ll we will find them, making the substi- 
tution fco = ifc4 in ([2]). 

In view of these two facts, the EM form factor can 
be expressed through the Euclidean BS amplitude only 
approximately. Below, we will study the accuracy of the 
following two approximations. 

(i) Naive Euclidean form factor. In this case, the form 
factor is found by the naive substitution fco — ikj^ in the 
Minkowski expression ([8]) . This corresponds to an approx- 
imate Wick rotation which disregards singularities. How- 
ever, the BS amplitude in the complex plane <?£;(fc4,fc') 
can be found exactly, by substituting in eq. ([2]) the com- 
plex values (fT5|l of boosted arguments. 

(m) Naive Euclidean form factor in the static approxi- 
mation. In this case, the form factor is still found by the 
substitution fco — ik^ in ([8]). In addition, the boosted 
amplitude <?^°**(fc4, fe;p) = (I>E{k'^,k') is approximately 
replaced by the amplitude at rest <P^£°'^*{k4,k; M,p — 
0) = <pE(k4,k). Due to that, the form factor is expressed 
through the Euclidean BS amplitude with real arguments. 



4.1 Naive Euclidean form factor 

In order to obtain the naive Euclidean form factor, we 
start with the Minkowski space formula ([8]). We use the 
Breit frame defined as: 



P 



P, Po = Po 



-fcn 



and shift the integration variable: fco — > —fco + ^Po- The 
spatial components of eq. ([5]) in the Breit frame are triv- 
ially satisfied (0=0). Taking the time-component, we get: 



2po^M(g' 



id^k 



Po 
2 



k' 



X {po + 2ko)<pM ( fco, - k;p \ $m ffco, -^P 



k-p' 



(16) 
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We simply replace: fco = 1/04 with real k^, that is, we 
neglect the contributions of singularities crossed by the 
rotated contour, and we obtain: 



F] 



maive /r^2 



2po(27r) 
,2 , 7.2 



(17) 



where f?^°**(A;4, A;;p) is defined in (HI]). Substituting in 
r.h.s. of ([13]) the BS ampHtude from eq. ([2]), one gets: 



1 /"OO 

dz / £^7 
47r J_i Jo 



(18) 



5(7,2) 



■ n3 



[7 + — + fc| + fc2 — {ipok4 — p ■ k) z ~ ie\ 



After substituting in p?|) . the form factor 
pnaive f^Q2-^ is cxprcssed as: 



1 



/ dz' g{i,z') / d7 



00 /'OO 



X / dzg{'j,z) / dfc4 / d'^k f{k4,k,p), 
-1 Jo Jo 

(19) 



where 



f{h,k,p) 



1 



Re[U{ki,k,p)] 



U{ki,k,p) 
with 



+ rn^ + (fc4 + iipo)^] {2iki+po) 



D ^ A- p ■ k{l + z) + - {1 + 2z)p^ - ikipaz - ie, 

A = j + m^- -M^ + kl + k^ 
4 



(20) 



and po = + p^ . D',A' are obtained from D,^ by 

the replacements 7 — > 7', z — > z'. The normahzation factor 
j^naive jg ^^^^^ fo^^^j fj.^^ ^-j^g Condition Fg'^™'=(0) = 1. 

Like Minkowski space form factor (fT3|) . the form fac- 
tor ([TO)) is expressed through the function (7(7,2), satis- 
fying ([3|), but differs from ([T3|) by neglecting singularities 
crossed when performing the Wick rotation. Their com- 
parison in sect. [Sj will show the error induced by this ap- 
proximation. 

Note that the denominators D and D', for some mo- 
mentum transfer and the integration variables, may be 
zero. For example, for 7 = 0, 2; = 1, A:4 = we get: 



(we used that Q ~ 4p ). Since m — |M is positive, D 
is always positive too if 

0^ < 4(4^2 _ j^.j2y (^21) 

If > 4(4to^ — M^), D crosses zero for some particular 
values of p and k. This singularity is, of course, integrable 
(the form factor is always finite), but it is a source of 
numerical instability. 



4.2 Naive Euclidean form factor in the static 
approximation 

As explained in the previous section, the form factor ([17p 
is expressed through the BS amphtude ^^°^*(A;4, k,p) which 
for p ^ 0, is in its turn expressed through the Euclidean 
BS amplitude in the complex plane. If we replace the 
latter by the Euclidean BS amplitude in the rest frame 
<pEiki,k) = ^^°''*(fc4, fc; M,p = 0), satisfying equation 
([4]), we obtain the form factor in the so called static ap- 
proximation, which reads: 



1 f dkjid'^k 



k' 



1 2 
4P0 



X [ki, \p - 'Pe {k4,-\p - k^ (22) 

Notice that no approximation in the kinematical factor is 
done, though we omit the odd degrees of k4, since after 
integration over dk^ they give zero. 

As mentioned, 'pE{k4,k) can be found from equation 
([1]). If 3(7, z) is known, <?£;(fc4, k) can be alternatively ob- 
tained by equation ([T8]) at p = 0. In this case, the integral 
(fTSl) is reduced to: 



^ n ,^ ■ [\ r , 9{l,z)2A{A^-iklM^z^) 
Jo Jo V47r(A2 -I- kjM^z^)^ 

where A is defined in ([^U)) . 



5 EM form factor via light-front wave 
function 

Knowing the Minkowski BS amplitude, we can find the 
light-front wave function [17j : 

^(fe^,^)^ (^^il^lif^ r^$,,(k + Pu;,p)df3. (23) 
t:{uj-p) J_^ 

Here w is a four-vector with = 0, determining the ori- 
entation of the light-front plane. The perp-components of 
vectors, which appear below, are defined relative to the 
direction u). Relation (|23p is independent of any model. 
Substituting ([2]) into ([23]) . we find the two-body light- front 
wave function: 



ip{k±,x) 



1 



D = m^- ^M^ - + (p - fc) 

4 lb 



x(l — a;)5(7, 1 — 2x)d^ 



(24) 
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The form factor is expressed through this wave function 
as (see e.g. [I7|): 



If d k 



<Pk 1 dx 



2x(l - xY 
(25) 

where Q\ = Q^. Substituting in the wave function 
■!/'(fe_L, x) determined by eq. (j24p and using the formula 

1 6w(l — u)du 



{au + b{l-u)^ 



we can easily integrate over A;_l and write the form factor 
as: 



Flfd{Q^) = 



1 



'CXD pOO pi pi 

/ / d-f dx du 

2^tt^Nlfd Ja Jo Jo Jo 



x(l - x) u(l - u) 5(7, 2x - l).g(7', 2x - 1) 
[m7 + (1 - m)7' + m(1 - u)x^Q^ +m? - x{l - x)M'^\ 



(26) 



naive Euclidean form factor calculated with boosted Eu- 
clidean BS amplitude by eq. Dashed curve denotes 
the form factor in the static approximation, eq. (j22p . The 
difference between solid and dotted curves shows that in- 
deed some singularities are missed and, therefore, the Wick 
rotation in the variable ko results in an inaccuracy. The 
static approximation (dashed curve) generates an addi- 
tional error, increasing with . 

Binding energy B = \ corresponds to M = m. In this 
case, the condition ([2T|l is violated if > Ylw?. Indeed, 
our numerical calculation became unstable if crosses 
K, Vlw?. That is why the domain of in fig. [5] does 
not exceeds Vdw?. This difficulty is absent in the static 
approximation. 

Figure [6] shows the comparison between form factors 
calculated in Minkowski space (solid) and in the static 
approximation (dashed) in a wider domain of momentum 
transfer. For high Q^, the static form factor is smaller 
than the Minkowski one by at least a factor 10. A zero 
in the static form factor at ~ 38 is an artefact of 
the static approximation, since it is absent in the exact 
(Minkowski space) form factor. 



6 Numerical results 

All the calculations given below have been done with the 
BS amplitude found for the ladder-l-cross ladder kernel. 
The constituent mass m = 1, the exchange mass = 0.5 
and the coupling constant a has been adjusted to provide 
the binding energy B — 1. 




(Q/m)' 




1,.05 

20 40 60 80 100 

(Q/m)' 

Fig. 6. Form factor via Minkowski BS amplitude (solid curve) 
and in static approximation (dashed). 

Figure [7] shows the comparison between the form fac- 
tors calculated in Minkowski space (solid curve, the same 
as in fig. ^ and in LED, eq. (dot-dashed curve). 

The LED form factor is almost indistinguishable from the 
Minkowski one in all domain of momentum transfer. 



Fig. 5. Solid curve: form factor via Minkowski BS amplitude. 
Dotted curve: form factor, after Wick rotation, with boosted — ^ , . 
Euclidean BS amplitude. Dashed curve: form factor in static ' »-OnCIUSIOn 
approximation . 

We have applied the solution of the BS equation in Min- 
kowski space, found by the method developed in [9, lOj, to 
For < 10 (" JLab domain"), the naive Euclidean calculate the EM form factor and to evaluate the inaccu- 
form factor and its static approximation are compared racy of different approximations available in the literature, 
with the Minkowski one in fig.O Solid curve is the Minkowski This method gives the BS amplitudes both in Minkowski 
space calculation, eq. (|13p . Dotted curve represents the and Euclidean spaces as well as in the full complex plane. 
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omagnetic form factor via Bethe-Salpeter amplitude. 




Fig. 7. Form factor via Minkowski BS amplitude (solid curve, 
the same as in fig. |6]) and in LFD (dot-dashed curve). 



We presented two additional validity tests of this method 
and demonstrated that it gives the same Euclidean BS am- 
plitude as the one found by directly solving equation 
For the ladder and ladder +cross ladder kernel, it gives 
the same binding energy. 

We calculated the electromagnetic form factor exactly, 
via Minkowski space BS amplitude. To express it through 
the Euclidean solution, one should carry out the Wick 
rotation, which, however, requires to incorporate the con- 
tributions of the singularities, crossed by the rotating inte- 
gration contour. In the naive Euclidean form factor, they 
are omitted. In addition, after Wick rotation, the Eu- 
clidean BS amplitude in a moving reference frame (i.e., 
boosted BS amplitude) is expressed through the rest frame 
one, depending on complex arguments. By our method 0, 
we find the BS amplitude in complex plane and analyze 
the error resulting from naive Wick rotation. The error 
increases with momentum transfer and at « 10 
(" JLab domain") it is about 30%. In the static approxima- 
tion, the error becomes larger, so that at 10 the 
Minkowski- and static-approximation form factors differ 
by one order of magnitude. The three form factors - the 
exact one from Minkowski BS amplitude, the Euclidean 
boosted one and in static approximation - are found to 
be close to each other (within a few per cents) only at 
relatively small momentum transfer < m? . 

The form factors calculated using the Minkowski space 
BS amplitude and the light-front wave function coincide 
with each other with very high accuracy. They are almost 
indistinguishable . 

Note that, in contrast to the Minkowski space BS am- 
plitude, the light-front wave function is not singular and 
can be found directly, from the corresponding 3D equation 
[T7], without using any BS formalism and eq. This 
advantage, together with accurate result for the form fac- 
tor (see figure [7]) , is one of the reasons which makes the 
application of the light-front approach to the EM form 
factor rather attractive. 



The system of spinless particles considered in this work, 
provides a simple model giving a lower limit of different 
approximations accuracy to the form factor. One can ex- 
pect that incorporating spin, these errors would increase. 
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